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ABSTRACT 

The  equations  governing  the  flow  of  viscoelastic  liquids  are  classified 
according  to  the  symbol  of  their  differential  operators.  Propagation  of 
singularities  is  discussed  and  conditions  for  a  change  of  type  are 
investigated.  The  vorticity  equation  for  steady  flow  can  change  type  when  a 
critical  condition  involving  speed  and  stresses  is  satisfied.  This  leads  to  a 
partitioning  of  the  field  of  flow  into  subcritical  and  supercritical  regions, 
as  in  the  problem  of  transonic  flow. 
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SIGNIFICANCE  AND  EXPLANATION 


Viscoelastic  liquids  exhibit  a  number  of  strange  and  unexpected  phenomena 
when  the  influence  of  elasticity  (the  "Weissenberg  number")  is  high. 

Numerical  simulations  of  high  Weissenberg  number  flows  have  encountered 
considerable  difficulties. 

This  paper  advances  the  idea  that  some  properties  of  flow  as  well  as  some 
problems  of  numerical  integration  may  be  associated  with  a  change  of  type  in 
the  governing  equations.  Many  models  of  viscoelastic  fluids  predict  such  a 
change  of  type  for  steady  flow,  similar  to  a  sonic  transition  in  gas 
dynamics.  We  investigate  the  conditions  for  such  a  transition  and  discuss 
various  applications. 
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phase  changes  and  may  be  relevant  In  analyzing  some  instabilities  in  the  extrusion  of 
polymers  from  capillary  tubes.  Chapter  4  discusses  characteristics  and  classification  of 
type  for  first  order  quasilinear  systems. 

In  Chapter  5,  we  look  at  constitutive  equations  for  viscoelastic  fluids  from  the  point 
of  view  of  classification  of  type.  For  this,  we  have  to  maintain  a  distinction  between 
fluids  with  and  without  Newtonian  viscosity.  In  Oldroyd  models,  the  term  with  Newtonian 
viscosity  is  the  one  associated  with  a  retardation  time.  The  addition  of  even  small 
amounts  of  Newtonian  viscosity  can  smooth  discontinuities,  replacing  sharp  fronts  by  thin 
layers  and  thus  masking  the  underlying  dynamics.  To  highlight  the  effect  of  hyperbolicity. 


* 


II 


we  confine  our  attention  to  models  without  Newtonian  viscosity.  In  particular,  we  focus  on 

a  three-parameter  family  of  nonlinear  Oldroyd  models  containing  the  upper  and  lower 

convected  and  corotational  Maxwell  models.  The  occurrence  of  Hadamard  type  instabilities 
for  these  models  is  discussed.  Theee  models  also  form  the  basis  for  the  discussion  of 
steady  flows  in  Chapters  9-11.  We  also  discuss  more  general  models  of  integral  type.  It 
is  shown  that  the  principal  part  of  the  linearization  at  any  given  motion  has  the  form  of  a 

rate  equation  not  involving  integrals,  provided  that  the  integral  kernels  have  sufficient 

smoothness.  Thus  the  discussion  of  change  of  type  does  not  necessarily  require  a  special 
constitutive  model. 

Chapter  6  discusses  the  linear  system  of  equations  for  motion  perturbing  rest.  The 
wave  speed  along  characteristics  is  given  by  /g(0)/p,  where  G(0)  is  the  instantaneous 
value  of  the  relaxation  modulus  G{s)  and  p  is  the  density.  We  review  recent  results  on 
the  propagation  of  slip  surfaces  for  velocity  and  displacement,  which  show  in  particular 
the  crucial  dependence  on  the  nature  of  the  kernel  G(s).  In  particular,  consideration  is 
given  to  the  possibility  that  G(0)  or  G'(0)  may  be  infinite.  In  Chapter  7,  we  discuss 
the  formation  and  propagation  of  slip  surfaces  in  nonlinear  shearing  problems  treated  by 
Coleman  and  Gurtin  [6],  [7]  and  Slemrod  [431,  [44].  We  discuss  the  application  of  their 
results  to  melt  fracture. 

In  Chapter  8  we  take  up  the  analysis  of  change  of  type  in  steady  problems.  This  is  a 
natural  question  from  a  mathematical  point  of  view,  but  the  first  studies  of  it  in  the 
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theory  of  viscoelastic  fluids  seem  to  he  in  the  work  of  Butkevich  [40,  p.  44-45],  who 
analysed  the  two  dimensional  equations  for  an  upper  convected  Maxwell  model.  Ultraan  and 
Penn  [49]  and  Luskin  [27]  classified  the  linearised  equations  perturbing  uniform  flow  with 
velocity  U  of  an  upper  convected  Maxwell  fluid.  Our  analysis  in  Chapter  8  generalises 
the  results  of  Ultman  and  Denn  and  Luskin  to  a  wider  class  of  constitutive  laws.  There  is 
a  change  of  type  leading  to  real  characteristics  when  the  viscoelastic  Mach  number 

M  -  U/c,  c  -  /G ( 0 ) /p 

exceeds  one.  Wie  vorticity  is  identified  as  the  variable  which  can  become  discontinuous 
along  these  characteristics.  We  shall,  somewhat  loosely,  say  that  "the  vorticity  changes 
type”.  In  Chapter  9,  we  give  a  complete  classification  of  the  quasilinear  system 
describing  the  upper  convected  Maxwell  model  in  arbitrary  steady  two-dimensional  motions. 
The  streamlines  are  double  characteristics.  The  vorticity  changes  its  type  when  the  speeds 
are  great  enough.  In  the  supercritical  (hyperbolic)  case,  there  are  two  families  of  real 
characteristics  for  the  vorticity,  but  the  formula  for  the  characteristics  depends  on  the 
solution.  There  are  also  complex  roots  to  the  characteristic  equation  associated  with  the 
elliptic  equation  giving  the  vorticity  as  the  Laplacian  of  the  stream  function.  In  Chapter 
10,  we  discuss  a  number  of  specific  flows  for  an  upper  convected  Maxwell  fluid.  These 
flows  include  plane  parallel  shear  flow,  steady  extensional  flow,  sink  flow  in  the  plane 
and  shear  flow  outside  a  rotating  cylinder.  We  discuss  characteristics  for  motions 
perturbing  those  flows  and  characterize  the  regions  of  flow  where  the  vorticity  equation  is 
hyperbolic.  In  Chapter  11,  we  extend  our  results  to  a  three  parameter  family  of  Oldroyd 
models  which  contains  the  upper  and  lower  convected  and  corotational  Maxwell  models  as 
special  cases.  The  vorticity  is  again  identified  as  the  variable  which  changes  its  type. 

We  compute  the  characteristic  directions  for  the  nonlinear  problem  without  approximation. 

We  exhibit  special  cases  which  show  that  the  partitioning  of  the  flow  into  sub-  and 
supercritical  regions  is  model  sensitive.  It  is  therefore  desirable  to  develop  this  type 
of  theory  on  a  high  level  of  generality,  suppressing  models.  We  take  some  steps  in  that 
direction  in  Chapter  12,  where  we  study  fading  memory  fluids  of  Coleman-Noll  type. 
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2.  THE  HIGH  WEISSENBERG  HUMBER  PROBLEM 


The  Weissenberg  number  is  e  number  which  is  large  when  the  elasticity  of  the  fluid  is 
important.  It  is  defined  differently  by  different  authors.  Many  authors  have  it 
proportional  to  a  relaxation  time  for  the  fluid,  or  it  may  be  viewed  as  the  ratio  of  the 
first  normal  stress  to  the  shear  stress.  The  parameter  of  interest  in  this  study  is  a 
viscoelastic  Mach  number,  which  is  the  ratio  of  a  typical  speed  to  a  wave  speed,  defined 
differently  for  each  problem,  and  is  large  when  the  relaxation  time  is  large. 

The  high  Weissenberg  number  problem  is  that  the  numerical  integration  of  equations  for 
viscoelastic  fluids  cannot  be  done  when  the  Weissenberg  number  is  high,  or  even  moderate. 
This  problem  is  not  solved  by  changing  the  model  or  by  different  methods  of  numerical 
integration.  Maybe  there  is  an  underlying  mathematical  problem. 

It  is  probable  that  some  numerical  problems  at  high  Weissenberg  number  are  actually 
associated  with  a  change  of  type,  like  the  transition  from  subsonic  flow  to  supersonic 
flow,  and  that  the  solution  of  the  problem  is  to  be  sought  in  various  hyperbolic 
algorithms,  especially  those  recently  introduced  for  transonic  flow. 

TO  compute  subsonic  flow  you  use  some  central  differences.  To  compute  supersonic  flow 
you  use  the  method  of  characteristics.  It  would  be  a  disaster  to  try  to  do  supersonic  flow 
by  central  differencing  of  the  type  used  for  taplaces  equation. 

In  the  flow  over  a  bump,  say  an  airfoil  with  the  free  stream  ever  so  slightly  less 
than  M  -  1 , 

.75  <M  <  1  , 

we  get  a  supersonic  bubble  with  unknown  boundaries  (see  Fig.  1). 


Fig.  1 
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To  solve  this  problem  you  have  to  find  the  sonic  line,  the  position  and  strength  of  the 
shock  wave.  This  Is  a  very  hard  free  boundary  problem.  It  wasn't  solved  until  1971  when 
Murman  and  Cole  [32]  realized  that  upwind  differencing  was  necessary  In  the  supersonic  part 
of  the  flow.  In  central  differencing  the  nodal  point  is  at  the  center.  In  upwind 
differencing,  the  information  at  a  nodal  point  is  determined  only  by  the  flow  upstream  (see 
Fig.  2). 

O  e 

0  *  0  0  0  * 

O  o 

central  differencing  upwind  differencing 

Fig.  2. 

Murman  and  Cole  studied  the  small  disturbance  equations.  They  derived  a  switching  scheme 
of  numerical  analysis  which  tells  the  computer  to  use  central  differencing  if  the  flow  is 
subsonic  and  upwind  differencing  if  the  flow  is  supersonic.  Their  upwind  differencing 
equation  can  be  interpreted  as  approximating  a  differential  equation  with  an  artificial 
viscosity  proportional  to  the  mesh  size  [31]. 

Murman  and  Coles'  method  was  the  first  success.  But  this  method  is  too  simple  for  the 
full  nonlinear  potential.  This  more  complicated  problem  was  successfully  attacked  by  the 
artificial  viscosity  method  of  Jameson  [19],  [20],  whose  work  makes  transonic  computation 
possible  in  a  practical  sense. 

People  doing  flow  computations  for  viscoelastic  fluids  are  also  able  to  go  to  higher 
Welssenberg  numbers  when  they  have  constitutive  equations  with  more  Newtonian  viscosity. 
This  procedure  masks  the  problem  of  dealing  with  change  of  type  instead  of  solving  it. 


3.  CONCEPTS  AND  SOME  APPLICATIONS  OF  CHANGE  OF  TYPE 


This  paper  deals  with  equations  which  undergo  a  change  of  type.  To  make  our  notions 
precise  we  shall  need  some  classical  definitions  related  to  the  type  of  a  partial 
differential  equation. 

Consider  the  linear  differential  operator 


(3.1) 

where 


p(x,t. 


3  3  3  n 

3t'  3x, . 3x  j 

1  n 


at  *  (X1  ,X2, . . .  ,xn)  and  t 
are  space  and  time  coordinates.  We  define  the 

(3.2)  Symbol  of  P  -  P<x,t,i50,ie., , . . .  ,i5n>  , 

_  3  3  3 

where  i  -  /-I.  To  form  the  symbol  we  replace  the  arguments  — ,  r — #  •  •  •  »r —  of  P 

d  t  dx4  d  x 

1  n 

with  the  Fourier  variables  i5g , iE^ , . . . ,i?n.  In  this  way  we  obtain  a  polynomial  in  the 

n  2 

real  variables  E.  Hie  symbol  of  the  Laplace  operator  -6  is  \  E,;  the  symbol  of  the 

32  2  n  2  1-1  3 

wave  operator  — r  -  A  is  “E.  +  5,  E.;  the  symbol  of  the  heat  operator  r—  -  A  is 

3t  i-1 

n  2 

i£  +  \  E  .  The  symbol  for  a  system  of  equations  is  defined  in  a  similar  fashion  and  is 

i-1  1 

a  matrix  with  polynomial  entries. 

Characteristic  curves  are  lines  along  which  discontinuous  data  may  propagate.  In 
dimensions  higher  than  two  we  may  speak  of  characteristic  surfaces.  Let  m  be  the  highest 
order  of  the  derivatives  in  P.  Then 


P  -  l  a  (x,t)3a  +  l  a  (x,t)3a 
|a|-m  |a|<m 


where  a  -  (Og.a^ 


The  equation 


, u^)  is  a  multi-index,  |a|  -  \  and 


,l«l 


,  a°a  a1  ,  “n 
3 1  3x  .  . • .  3x 

1  n 


l  a  (x,t)aa  =0,  a  =  (o  ,...,o  )  , 


(3.3) 


■SS 


w 


ia  called  tha  characteristic  equation  for  P,  Only  the  principal  part  of  P,  the  terms  of 


highest  order,  appears  in  (3.3). 


A  surface  S  in  (jj,t)  space  is  characteristic  for  P  at  a  point  s  e  S  if  the 


normal  vector  to  S  at  s  satisfies  the  characteristic  equation.  If  a  =  (a. 


,cr  )  is 


a  unit  normal  vector  at  s,  S  is  characteristic  for  P  if  and  only  if 


X  =  1  and  \  a  (x)oa  -  0. 

0  |a|“m 


n  2  2  n  2 

The  characteristic  equation  for  Laplaces  equation  \  3  u/3x  “  0  is  \  o  = 


There  are  no  real  characteristics  because  (3.4).  is  not  satisfied.  More  generally,  t 


operators  P  for  which,  at  every  point  (j{,t)  the  equation  (3.4)2  has  no  nontrivial  ce 


zeroes  are  called  elliptic.  For  systems,  ellipticity  means  that  the  only  real  zeroes  of 


the  determinant  of  the  matrix  symbol  A(x,€ 1 , • . . ,5  )  are  (E ^ ,S2 / • •• «£„)  =  (°»c 


Elliptic  problems  have  existence,  uniqueness  and  continuous  dependence  on  data  (are 


well  posed)  as  boundary  value  problems  [1],  [26]. 


The  initial  value  problem,  the  Cauchy  problem,  is  not  well  posed  for  elliptic 


equations.  For  example,  a  Cauchy  problem  for  Laplaces  equation  in  the  domain 


D  »  {x,y»  x  >  0,  -•  <  y  <  «•}  is 


Au  *  0  in  D  , 


u(0,y)  -  0  , 


fjj  (0,y)  -  U(y)  , 


U(y)  “  —  sin  ny,  p  >  0  . 


The  solution  of  (3.5)  is 


u(x,y)  “  sin  ny  sinh  nx  . 

n 


The  mapping  (u,  ■f^)jx.o  — u  for  x  >  0  is  not  continuous  since  U(y)  is  small 


n  is  large  and  u(x,y)  is  very  big.  Small  data  at  x  ■  0  lead  to  larger  and  larger 


<•.  • .  - 
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oscillations  for  x  >  0.  This  lack  of  continuous  dependence  is  called  Hadamard 
Instability.  It  can  be  shown  that  (3.5)  has  no  solution  if  U( * )  is  not  analytic. 

The  initial  value  problem,  or  mixed  initial-boundary  value  problems  are  well-posed  for 
hyperbolic  equations  like  the  wave  equation.  For  example,  the  characteristic  equation 
( 3.4)j 


a 


2 

0 


0 


for  the  n-dimensional  wave  equation 


9  u  2. 
- -  c  Au 

at2 


//  c2  +  1. 


satisfies  the  characteristic  equation  (3.4)1  when  <JQ  -  ±c //  c“  +  1 .  Therefore  a  surface 
is  characteristic  for  the  wave  equation  if  and  only  if  its  normal  makes  an  angle 

8,  cos  8  »  c //  c2  +  1 , 

2  ,.„2 


with  the  t  axis.  For  the  one-dimensional  wave  equation 

A  “  3*/3xz,  this  implies  that  the  family  of  lines  x±ct  «  const  are  characteristic. 

The  operator  P  of  (3.1)  is  called  strictly  hyperbolic  if  all  the  roots  of  the 

principal  part  of  its  symbol  (3.2)  are  real  and  distinct  for  all  (£,,...,£  )  e  Vp\0.  The 

»  n 

Cauchy  problem  is  well  posed  and  the  boundary  value  problem  is  ill  posed  for  hyperbolic 
equations.  The  backward  Cauchy  problem  where  t  is  replaced  with  -t  is  also  well-posed 
for  hyperbolic  equations. 

The  Cauchy  problem  is  well  posed  for  parabolic  problems  but  the  backward  Cauchy 
problem  is  ill  posed.  The  classic  example  of  a  parabolic  equation  is  the  heat  equation, 
3u/3t  =  Au.  The  characteristic  equation  (3.4)2  is 

n  ? 

l  O2=0. 

?.«1 

2 

Hence,  from  (3.4)^,  ”  1  and  the  characteristic  surfaces  are  the  hyperplanes 

t  •  const.  The  Cauchy  problem  is  not  well  posed  for  the  backward  heat  equation  ^  «  -Au. 

<7  t 

Operators  of  the  form  +  Lu,  where  L,  like  -A,  is  a  positive  definite  elliptic 

o  t 


operator,  are  parabolic.  These  operators  are  strongly  dissipative  and  lead  to  diffusion 


rather  than  to  propagation.  Unlike  hyperbolic  operators,  parabolic  operators  will  smooth 
initially  discontinuous  Cauchy  data. 


Two  homogeneous  scalar  operators  are  said  to  be  of  the  same  type,  if  up  to  a 
transformation  of  the  independent  variables,  their  symbols  have  the  same  asymptotic 
behavior  at  infinity.  If  the  asymptotic  behavior  of  the  symbol  changes,  then  we  say  that 
the  equation  changes  type.  For  example,  the  Tricomi  equation 


is  hyperbolic  when  y  <  0 
system 


32u  32u 

y  — 2  — 2 

Sx*  Sy*4 


and  elliptic  when  y  >  0. 


0 


Another  example  is  the  quasilinear 


3u  3<r(v) 

3t  “  3x  ' 

(3.6) 

3v  3u 
3t  ”  "Sx' 

which  is  hyperbolic  for  o' (v)  >  0  and  elliptic  for  o'(v)  <  0.  These  problems  all 
involve  a  change  in  the  sign  of  the  symbol  and  Hadamard  instabilities.  If  we  start  with 
initial  data  in  the  hyperbolic  region,  and  if  the  solution  of  the  Cauchy  problem  enters  the 
elliptic  region,  one  has  Hadamard  instability  and  a  loss  of  evolutionarity . 

Problems  of  the  form  (3.6)  suggest  models  for  theories  of  phase  changes  in  solids  and 
fluids.  The  van  der  Waals  gas  is  a  well-known  classical  example.  In  solid  mechanics, 
ideas  of  this  type  were  introduced  by  J.  Ericksen  [12]  in  his  study  of  elastic  bars.  We 
may  suppose  that  the  graph  of  o(v)  is  as  shown  in  Fig.  3. 


32v  3  ,  , 

—  “  37 


3t 


(v)  gn 


The  system  (3.6)  is  hyperbolic  when  o’(v)  >  0. 

The  elliptic  branch  is  unstable  in  the  sense  of  Hadamard 


The  solid  lines,  where  c'(v)  >  0,  lead  to  a  hyperbolic  equation  and  the  dashed  line  leads 
to  an  elliptic  equation.  The  elliptic  portion  is  _ ejected  because  it  will  exhibit  Hadamard 
instabilities)  and  actual  solutions  are  required  to  operate  only  on  the  hyperbolic  parts  of 
the  curve.  •  This  leads  to  spatially  segregated  solutions,  separated  by  lines  of 
discontinuity,  each  part  operating  on  a  different  hyperbolic  branch  of  the  curve.  There  is 
hysteresis  and  abrupt  transitions  in  the  response  of  such  models.  These  features  are  all 
present  in  the  recent  study  of  Hunter  and  Slemrod  [17],  which  attempts  to  explain  some 
observations  of  Tordella  [47]  of  a  type  of  melt  fracture  called  ripple.  This  phenomenon 
shows  hysteresis  loops,  double-valued  shear  rates  at  certain  stresses  and  spatially 
segregated  flow  regimes.  Similar  ideas  have  also  been  used  to  explain  to  phenomenon  of 
necking  occurring  in  cold  drawing  of  polymers  [8]. 

Regirer  and  Rutkevich  [36]  have  considered  fluids  of  the  Reiner-Rivlin  type  which 
exhibit  change  of  type.  Their  constitutive  law  is 

t  -  -pi  +  nftiDD  , 

It  2 

where  D  ■*  —  (Vu  +  (7u)  ),  11  »  tr  D  .  Written  in  terms  of  a  stream  function 

U  «  (u,v)  «  (♦  ,-♦  ),  the  equation  governing  steady  two-dimensional  flows  is  as  follows, 

y  ^ 


(3.7) 


def 
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where  H(t|>)  is  a  nonlinear  thrid  order  operator  and  the  coefficients  a  are  nonlinear 
functions  of  the  second  derivatives  of  The  characteristic  curves  y(x)  are  solutions 

of 

4  3  2 

(3.8)  a.y  +  4a, y  +  2a, y  -  4a, y  +  a,  *=  0  . 

lx  3  x  2  x  3  x  1 

It  turns  out  that  there  are  the  following  three  cases 

(i)  f  +  211  f'  >  0  (no  real  roots,  elliptic), 

(ii)  f  +  2X1  f'  -  0  (parabolic), 

(iii)  f  +  211  f'  <  0  (four  real  roots,  hyperbolic). 


The  hyperbolic  regions  are  those  where  the  stress  decreases  as  a  function  of  shear  rate, 


and  the  elliptic  regions  are  those  where  it  increases.  The  unsteady  problem  corresponding 
to  (3.7)  is 

P  1^-  (&*)  -  Zi(l|/ )  -  H(t|>)  . 

When  the  right  side  is  elliptic,  this  problem  is  parabolic  and  evolutionary  (see  Gelfand 
(14]).  When  the  right  side  changes  type,  evolutionarity  is  lost  and  Hadamard  instability 
occurs.  Changes  of  type  and  Hadamard  instabilities  can  occur  in  rheological  problems  which 
are  not  esssentially  one-dimensional  and  need  not  necessarily  be  associated  with  non- 
monotone  shear  stress-shear  rate  laws.  An  interesting  case  of  this  type  arises  in  a 
stability  analysis  of  plane  Couette  flow  by  Akbay,  Becker,  Krozer  and  Sponagel  [3],  which 
they  suggest  as  an  explanation  for  melt  fracture.  In  order  to  obtain  a  manageable 
equation,  they  introduce  the  "short  memory  approximation".  This  means  that,  in  the  memory 
integrals  occurring  in  the  equation  for  the  disturbances,  only  terms  of  first  order  in  the 
relaxation  time  of  the  fluid  are  kept.  Proceeding  thus,  they  find  the  following  linearized 
equation  for  the  stream  function  in  two  dimensions > 


p(k +  kx2  •  (n;  -  7l) 


-?1  .!■»  +TT.2»  »ii  . 
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m  m  0 


at  x2  »  0,  x2  •  h  . 
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Here  L  denotes  the  operator  — —  -  — — .  The  problem  is  posed  in  the  strip  -“  <  x.  <  00 , 

3x‘  3x‘ 

0  <  x2  <  h.  k  is  the  shear  rate  of  the  basic  Couette  flow  and  t(k),  N1 (k )  are  the  shear 
stress  and  the  first  normal  stress  difference  as  functions  of  the  shear  rate.  Akbay  et  al. 


find  that  (3.9)  admits  exponentially  growing  solutions  if 


K 


It  was  pointed  out  by  Ahrens,  Joseph,  Renardy  and  Renardy  [2]  that  this  instability  is 
associated  with  a  change  of  type.  If  we  consider  the  symbol  of  the  differential  operator, 
i.e.  if  we  formally  set  tt  "  o,  ”  ia,  — -  i8»  then  the  left  side  of  (3.9)  becomes 

u  L  O  Xi  dXe 


(3.11) 


p (a  +  KXjia) ( -a  -  B  ) 


and  the  right  side  becomes 


(3.12) 


-(n;  -  ^)<a2  -  e2)aB  +  t  '  (a2  -  62)2  +  ~  a2fi2. 


This  homogeneous  fourth  degree  polynomial  is  positive  definite  for  Me  <  4,  but  indefinite 
for  We  >  4.  For  We  >  4,  one  thus  expects  short-wave  instabilities  of  a  catastrophic 
nature,  i.e.  Re  c  becomes  arbitrarily  large  as  the  wave  length  tends  to  zero.  A  more 
precise,  rigorous  analysis  and  some  criticisms  of  the  assumption  of  short  memory  are  given 
in  (21. 
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4.  gOASXLXNEAR  SYSTEMS 

The  analysis  of  the  equations  of  viscoelastic  flow  will  be  framed  in  terms  of  systems 
of  equations  of  first  order.  We  consider  linear  systems  and  quasilinear  systems.  We  write 
the  quasilinear  system  as 


n  3u 

J*  St  iT"  “  S'  5  "  (t,xvx2,  ...,xn) 
1-0  1 


where  u  -  (uj,u2, . . .  ,Ujc)  is  a  k  vector  and  are  k  x  k  matrices  which  like  f  may 

depend  on  x^  and  on  the  components  of  u.  If  is  independent  of  u,  and  f  «  B*  u , 

then  (3.7)  is  a  linear  system. 

The  following  definitions  apply  to  both  linear  and  quasilinear  systems.  A  surface 
S  defined  by  the  equation  ♦(t,x1 , . . .  *xn)  -  is  characteristic  with  respect  to  (4.1)  at 


5  "  (t'x1 . V  if 


d«t(  l  At  |£-)(x)  -  0  . 


1-0  t 


If  ♦  -  xn  -  f(x0,...,xn_1),  then 

,4.3)  det(An  -  T  *,  H;)  *  °  • 

1-0  1 

Anyone  of  the  n  +  1  quantities  3$/9xj  in  (4.2)  may  be  regarded  as  an  eigenvalue,  we 
shall  say  (4.1)  is  hyperbolic  provided  that  4  -  is  non-singular  and  for  any  choice  of 
the  real  parameters  (Xf,  1  -  0,1, ,..,ni  1  *  u),  the  roots  a  -  a.  of 
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(4.4)  det(oA  -  l  X.4.)  -  0 
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t*U 

are  real  and  are  associated  with  k  linearly  independent  characteristic  vectors  v 


satisfying 
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First  order  »y» terns  can  be  of  mixed  type  with  real  and  complex  eigenvalues,  neither 
totally  elliptic  or  totally  hyperbolic. 

We  are  interested  in  two  dimensional  quasilinear  problems  of  the  form 

3u  3u  3u 

<4-6>  *  si  +  1  37  +  £  37  “  i  ■ 

We  consider  one-dimensional  evolution  problems  in  which  £  -  0  and  steady  problems  in 
which  £  “  0.  For  evolution  problems  we  suppose  that  A  is  not  singular  and  write 

3u  3u 

<4-7>  3?+8^-£‘ 

The  characteristic  surface  is  6(x,t)  -  0  and  (4.2)  becomes 

(4.8)  det(f±  1  ♦  B  |i)  -  0  . 

On  +(x,t)  >0,  we  have 

d*  -  |i  dx  +  |A  dt  -  °  . 

Hence,  (4.8)  may  be  written 

(4.9)  det(|  "  i)  "  0  < 
where  dx/dt  is  the  slope  of  the  characteristic. 

A  linear  system  of  the  form  (4.7)  is  said  to  be  of  evolution  type  if  |  has  only  real 
eigenvalues,  hyperbolic  if  of  evolution  type  and  B  can  be  made  diagonal,  strictly 
hyperbolic  if  B  has  simple  real  eigenvalues. 

If  the  Cauchy  problem  for  (4.7)  with  f  -  0  is  well-posed,  then  it  must  be  of 
evolution  type.  Solutions  of  the  form  u(x,t)  “  Be^*t+llx*  are  bounded  for  large  |t|  if 
and  only  if  the  eigenvalues  \/v  of  B  are  real  (Gelfand  [14]). 

Suppose  u(x,t)  is  given  on  a  curve  *(x,t)  «  0.  If  this  line  is  characteristic  then 
the  equation 

fa  .  ll  at  3a  dx 

ds  3t  ds  Tx  ds  ' 


where  t(s),  x(s)  is  a  parametric  representation  for  the  curve  $  «  0,  and  the 
quasilinear  equation  (4.7)  cannot  be  uniquely  Bolved  for  the  2k  derivatives  3u/3t  and 
3u/3x.  This  special  condition  requires  that  the  determinant  of  the  coefficients  of  the 


derivative*  vanieh 
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Identical  considerations  apply  for  the  quaailinear  steady  problem 

3u  3  u 

B  t—  +  C  t —  ■  f  . 

”  3x  »  3y 

Such  problems  are  frequently  associated  with  a  change  of  type,  like  transonic  flow,  in 
which  sous  regions  of  flow  are  subcritical  and  some  supercritical.  A  typical  example  is 
Tricomis  equation.  Other,  more  relevant  examples  are  derived  in  $8-11. 

It  is  not  always  possible  to  assign  a  definite  type  to  a  system  of  quaailinear 
equations.  He  may  have  both  real  and  complex  eigenvalues.  Nonlinear  problems  of  mixed 
type  have  not  been  thoroughly  studied  by  mathematicians.  Some  special  results  have  been 
given  by  Mock  [30] .  Here  it  is  perhaps  useful  to  give  some  simple  examples  from 
hydrodynamics . 

Consider  first  the  Euler  equation*  for  the  flow  of  inviscid,  incompressible  fluids  in 
two  diMnsions 

(4.10)  p(jj»V)u  +  Vp  «*  £,  div  jj  ■  0  . 

Let  (u,v)  be  the  components  of  u  with  respect  to  x  and  y.  Then  we  can  write  (4.10) 


(4.11) 


-  f  , 


where 


•r*K  +  &2-3y 


Tha  characteristic  equation  for  (4.11)  ia 


Hence 


dy  v 
dx  u 

±i 

dx 


.  streamlines  are  characteriatic. 


Tha  preaence  of  imaginary  roots  means  that  (4.11)  is  not  hyperbolic.  It  is  not  elliptic 
because  the  determinant  of  the  matrix  symbol  of  (4.11) 


u5t  +  ve2 


det 


u51  +  *2 


PC. 


-(uC,  +  vf-2)(f^  +  C2) 


vanishes  for  u£1  +  v?2  -  0. 

Another  example  is  from  the  theory  of  irrotational  water  waves.  In  this  case  the 
velocity  potential  is  elliptic  but  the  height  function  is  governed  by  a  hyperbolic  equation 
giving  rise  to  water  waves. 
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5.  CONSTITUTIVE  EQUATIONS 

A  constitutive  equation  relates  stress  to  deformation.  The  stress  in  viscoelastic 
fluids  depends  on  the  history  of  the  deformation.  Usually  the  history  is  defined  on  some 
strain  measure.  The  stress  in  Newtonian  fluids  depends  on  the  instantaneous  value  of  the 
velocity  gradient,  it  does  not  depend  on  the  history  of  the  deformation  but  only  on  the 
instantaneous  value  of  the  symmetric  part  of  the  velocity  gradient.  In  a  viscoelastic 
fluid  the  fluid  has  instantaneous  elasticity.  Elasticity  is  present  also  in  inviscid 
compressible  fluids.  For  unsteady  problems  elasticity  is  associated  with  hyperbolic, 
rather  than  parabolic  response.  In  the  elastic  case  discontinuities  propagate  along 
characteristics  which  do  not  exist  in  the  viscous  (Newtonian)  case.  It  is  necessary  to  be 
more  precise  about  the  difference  between  an  elastic  and  viscous  response. 

Many  constitutive  models  have  been  proposed.  Each  one  leads  to  different  answers  for 
the  same  problem  though  same  groups  have  similar  qualitative  properties.  In  problems  of 
changing  type  the  linearized  part  is  of  primary  importance.  The  linear  part  may  be  of 
three  types: 

1)  Constitutive  equations  with  some  viscosity.  The  viscosity  which  we  have  in  mind 
is  that  which  rheologists  sometimes  associate  with  a  retardation  time. 

2)  Constitutive  equation  without  "viscosity".  Integral  type  constitutive  equations 
with  smooth  kernels,  and  various  types  of  rate  equations  in  the  class  called  Maxwell  models 
are  of  this  type.  These  kind  of  equations  allow  propagation  of  rather  than  smoothing  of 
discontinuities.  In  some  nonlinear  models  [161,  [281,  [43],  [44],  [51]  discontinuities  may 
arise,  as  do  shocks  in  gas  dynamics,  from  smooth  data. 

3)  Integral  type  constitutive  equations  with  singular  kernels.  These  are  in  a  sense 
Intermediate  between  1)  and  2).  Depending  on  the  type  of  the  singularity,  the  wave  speeds 
may  be  finite  or  infinite.  However,  even  if  they  are  finite,  i.e.  real  characteristics 
exist,  there  is  no  propagation  of  discontinuities  (see  Chapter  6). 

The  stress  in  an  incompressible  fluid  is  given  by 
(5.1)  I  “  "Pj  +  l  , 

where  T,  the  determinate  stress  [50],  p.  176  (sometimes  called  the  extra  stress  [48]), 
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may  be  related  to  the  deformation,  whilst  p,  the  reaction  pressure,  is  determined  only 


through  the  equations  of  motion.  An  example  of  r  in  the  class  1)  of  constitutive 
equation  with  some  viscosity  is  the  Jeffreys  model  with  two  time  constants,  a  relaxation 
time  and  a  retardation  time  Xj.  This  model  may  be  written  in  rate  form 

3t  3a 

(5.2)  I  *  x,  -  hA  ♦  nx2  5=  . 

T 

where  n  is  a  constant,  called  the  zero  shear  rate  viscosity,  a  ■  7u  +  7u  ,  u  “  u(x,t), 
or  in  integral  form 

nX  x  t 

(5.3)  T  “  — r—  A  [u(  x,  t )  ]  +  —  ( 1  -  )  /  A[u(x,T]exp( - r — — )  dr  • 

The  constant  nX^/X^  is  a  second  viscosity  which  is  equal  to  the  zero  shear  rate  viscosity 
when  Xj  ■  X^.  It  is  this  viscosity  that  we  have  in  mind  when  we  talk  about 
"with"  or  "without"  viscosity.  In  $6  we  show  that  (5.2)  and  (5.3)  enjoy  a  certain  general 
status  when  they  are  regarded  as  holding  only  in  motions  which  perturb  a  state  of  rest. 

The  Maxwell  model  arises  from  (5.2)  and  (5.3)  when  X^  is  put  to  zero.  In  $6  we  note 
that  the  Maxwell  model  permits  propagation  of  waves  along  characteristic  with  a  finite 
velocity  of  propagation,  but  such  propagation  cannot  occur  for  Jeffreys  models;  more 
precisely  the  viscosity  term  X2n/X 1  smooths  discontinuities  in  the  same  way  that 
viscosity  smooths  the  discontinuities  of  solutions  of  Eulers  equations  for  an  lnviscid 
fluid. 

Nonlinear  models  can  be  classified  according  to  the  type  of  their  linearization. 
Popular  models  of  the  rate  type  include  those  due  to  Oldroyd  [34],  Leonov  [25],  Giesekus 
[15].  These  models  generalize  both  Maxwell  and  Jeffreys  type  fluids,  i.e.  some  have 
"viscosity",  some  do  not.  Popular  models  of  integral  type  include  K-BKZ  single  integral 
models  [4],  [22]  which  may  be  of  type  2)  or  3)  depending  on  the  nature  of  the  kernel,  and 
the  model  of  Curtiss  and  Bird  [11],  which  contains  a  viscosity  term. 


We  ahall  do  soma  work  with  a  three  constant  rate  equation  of  Oldroyd  type 


DT 

(5.4)  X  d£  +  i  "  n-  • 

Hare  O/Dt  is  an  Invariant  time  derivative 


(5.5) 


DT  3t 

—  -  +  (u'V)T  +  TO  -  Ox  -  a(DT  +  TD)  , 

Dt  St  *  -  “ 


where  -1  <  a  <  1,  g  “  j 

has  a  -  1  and 

(5.6) 

where  (Vu)^  ”  iu^/9*j. 


(7u  +  VuT),  2  “  J  ~  The  upper  convected  Maxwell  model 


DT  3x 

-  m  -  +  (u*V)T  -  VuT  “  TVu  i 

Dt  St  -  — 

The  lower  convected  Maxwell  model  has  a  ■  -1  and 


DT  3t 

(5.7)  5^  ”  +  (u*V)x  +  x7u  +  Vu  T  . 

The  corotational  Maxwell  model  has  a  “  0.  The  integral  model 


(5.8) 


t  , 

|  exp[-(t  -  T)/))(jJ'(T)  -  1]dx 


is  an  alternative  form  of  the  upper  convected  Maxwell  model  (a  »  1).  The  rate  form 
(5.4,6)  may  be  obtained  by  differentiating  (5.8)  partially  with  respect  to  t,  holding  x 
fixed.  The  expression 


(5.9) 


t 

I  exp(-(t  -  x  )/X]  [1  -  gt(x)]dx 


is  equivalent  to  (5.4,7)  with  a  “  ~1,  in  the  same  way. 

Rutkevich  (40],  [41]  studies  differential  constitutive  models  (5.4)  of  Oldroyd  type. 

He  linearises  these  equations  and  the  equations  of  motion  at  a  state  of  no  motion  (u  =  0) 
and  constant  stress.  He  finds  that  a  change  of  type  leading  to  imaginary  wave  speeds  and  a 
Hadamard  type  instability  occur  if  the  principal  values  of  x  satisfy  certain 
inequalities.  If  we  denote  these  principal  values  by  ti>t2>t3'  then  Butkevich's 
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instability  criterion  reads  as  follows: 


T3  <  "  T  '  *  “  1  » 


(5.10) 


a  -  -1, 


t.  -  t,  > 


1  3  T- 


In  order  to 


the  significance  of  these  results,  it  is  necessary  to  look  also  at 


the  integral  nwdels  (5.8)  and  (5.9)  corresponding  to  (5.4).  Since  and  ^  are 

positive  definite  matrices,  we  find  that  >  -  y  tor  a  “  1  and  T1  <  for  a  “  -1. 


Thus  Rutkevich’s  instability  criterion  can  not  be  achieved  and  must  be  considered 


unphysical.  This  was  also  noticed  by  Crochet  [10]  who  proved  directly  from  the 


differential  equation  that  T3  >  ~  ^  or,  resp.,  for  all  time,  if  this  is  the  case 


initially. 


Renardy  [39]  has  shown  that  in  fact  loss  of  evolutionarity  never  occurs  for  a  «  ±1 . 


These  two  models  are  special  cases  of  the  Kaye-BKZ  model  [22],  [4],  which  has  the  form 


(5.11) 


l  m  _[  a<t  “  St1(s)  "  3^  £tl8)^ds 


Here  a  is  a  positive  kernel  and  the  "strain  energy"  W  is  a  scalar  function  of 


I1  «  tr  (s)  and  I2  "  tr  ^(s).  Renardy  proved  that  the  initial  value  problem  is 


always  well  posed,  i.e.  Hadamard  type  instability  cannot  occur,  if  W  satisfies  a  strong 


ellipticity  condition  of  the  same  form  as  in  nonlinear  elasticity.  This  condition  is  in 


particular  satisfied  if  W  is  monotone  in  both  arguments  and  a  convex  function  of  /i^ 
and  /i^.  ttiis  obviously  includes  the  cases  W  «  Xj,  Vf  ■  I2,  corresponding  to  (5.8), 


(5.9).  No  such  objections  can  be  raised  in  the  corotational  case,  a  -  0.  in  fact,  we 


need  only  restrict  attention  to  motions  where  (J  -  0,  and  j  and  g  are  spatially 


homogeneous,  to  see  that  (5.4)  allows  -  x3  to  reach  any  value.  The  instability 


this  case  is  therefore  genuine. 
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It  is  instructive  to  look  at  special  flow  geometries.  For  time-dependent  simple 
shear-flow,  we  have  [21] 


)dt' )ds  . 


I  =  /  e  ^ ( s)  cos(/  1  -  a^  ]  >c  ( t  *  ) c 

s 

Differentiating  this  with  respect  to  t,  we  find 

t  «  tc(t)[y  -  /  ^  e  (t  S^//^lc<a)i/  1  -  a^  sin(/  1  -  a2  /  k  ( t'  )dt 1 )  ds] 


+  lower  order  terms  . 

A  change  of  type  occurs  when  the  expression  in  brackets  changes  sign.  For  |a|  <  1,  it  is 
possible  to  construct  histories  for  the  shear  rate  x(s)  such  that  this  is  the  case.  For 
steady  shear  flows,  x  =  kq,  the  expression  in  brackets  is  positive  and  there  is  no  change 
of  type,  even  though  the  shear  stress-shear  rate  law  for  a  *  ±1  is  non-monotone.  The 


case  of  simple  elongation  was  considered  by  Renardy  [38].  He  finds  that  a  change  of  type 

1 


can  occur  for  -1  <  a  <  -j.  Although  Rutkevich's  instability  for  a  “  ±1  is  unphysical  it 
is  nevertheless  relevant  for  numerical  calculations  (see  Crochet  [10]),  The  reason  is  that 


at  high  Weissenberg  numbers  the  eigenvalues  of  x  can  be  arbitrarily  close  to  the  stability 
n 


boundary  ±  j-  Numerical  errors  can  push  them  beyond  thiB  limit,  with  disastrous  results. 

Host  of  the  constitutive  equations  which  have  been  proposed,  all  the  ones  considered 
here,  are  simple  fluids  in  which  the  stress  is  determined  by  the  history  of  the  relative 
gradient  of  the  deformation 

£t<x,x)  “  vXt(*»T>  > 

where 

is  the  velocity  of  the  particle  x  •  jjt<x, t)  at  times  x  <  t.  The  constitutive  equation 
satisfying  material  objectivity  may  be  expressed  by  a  functional  [48, p.  80] 


(5. 12) 


x ( x, t )  -  F  C£t<2S'T >) 
T=-» 

on  the  history  of  the  relative  Cauchy  strain 


Ct(x,T)  -  FtFt  . 
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In  order  to  give  a  precise  meaning  to  an  expression  such  as  (5.12),  it  is  necessary  to 


specify  the  set  of  arguments  on  which  F  is  defined,  called  the  domain  of  F  or  dom  F. 

Saut  and  Joseph  [42],  extending  ideas  of  Coleman  and  Noll  [9],  have  proposed  to  classify 

constitutive  laws  according  to  the  choice  of  this  domain.  Roughly,  this  associates  the 

nature  of  the  stress-strain  relation  with  the  deformations  which  are  allowed.  In  the 

linearized  case  the  stresses  which  are  allowed  are  functionals  in  the  topological  dual  of 

dom  F.  The  smaller  dom  F  is,  the  larger  is  the  dual,  i.e.  the  more  constitutive  laws  are 

allowed.  To  be  more  precise,  let  some  given  history  and  g(x )  =  c(x,t  -  T)  be 

0 ,  , 

a  small  perturbation.  The  linearization  of  £  relative  to  C  (t  )  can  formally  be  written 

as  an  integral 

(5.13)  rjc^ft)  |  c  ( T )  ]  *  /  K(C°(T),t  -  T)c(T)dT  . 

—00  — 

T— * 

The  class  of  admissible  kernels  g  depends  on  which  deformations  c(t  )  are  allowed.  If, 

2 

following  Coleman  and  Noll,  g  is  restricted  to  a  weighted  L  -space,  then  g  must  also  be 
2 

in  a  weighted  L  -space.  If  g  is  restricted  to  a  Sobolev  space,  then  Dirac  measures  and 
derivatives  of  Dirac  measures  can  be  included  in  K,  thus  admitting  Jeffreys  type  models. 

To  see  the  effects  of  hyperbolicity  clearly  we  exclude  these  cases  and  adopt  (5.13) 
with  a  smooth  kernel  as  the  basis  for  our  study  of  change  of  type. 

It  will  be  useful  in  what  follows  to  specify  the  quantities  in  (5.13)  more 
precisely.  Let 

(5.14)  &t<x,T)  "  C°(x,T)  + 
be  the  particle  path  for 

x  =  £  (x,t)  =  5°(x,t)  , 

where  5°  is  the  position  of  x  in  some  given  motion  and  £(x,t)  is  a  perturbation  with 
E(x,t)  ■  0.  Quadratic  and  higher  order  terms  are  neglected.  To  this  order  of 
approximation 


C  (X,T)  =  C^(X,T)  +  C(X,T) 


(5.15) 


where 


(5.16)  C^(x,t)  “  (vi°)T75°f  C°(x,t)  -  1  , 

c(x,x)  -  <75°)TVC  +  (VC)TVe°,  c(  x , t )  -  0. 

Let  the  perturbed  extra  stress  be  denoted  by  T  and  the  unperturbed  extra  stress  by 
The  total  extra  stress  Is  +  r,  where 


(5.17) 


TQ  -  P  (c£(x,T)]  , 

T  -  /  K[C°(x,T),t  -  T]c(x,T)dT  . 


We  assume  that  £  is  a  smooth  function.  Applying  the  substantial  derivative 
+  u°*V),  where  u°  is  the  unperturbed  motion,  we  obtain 


(5.18) 


Roting  that  the 
zero,  i.e. 
(5.19) 


(|t  +  u°*7)*  “  S  +  u°*7)c(x,T)dT 

+  terms  of  lower  differential  order, 
position  at  time  t  of  a  particle  must  have  the  substantial  derivative 

<*Xt 

- - Jt—  +  «(x,t).VxXt(x,T)  -  0  , 


we  find 
(5.20) 


u°(x,t)*V  C° 


0  , 


2ft  0 

—  +  u  (x,t)»V  £  +  v(x,t)»V  £  «  0  - 

3t  x-  -  -  x- 

This  yields 

(ft  +  (u0<7))c  -  (|j  +  <u0»7))(V£0TV£  +  75TV5°) 

■  V5°T7{(|^  +  (u°»V))c}  +  (V{(|1  +  u°.v)5})TVC° 

+  lower  order  terms 


-7?ui.7(v*V5'')  -  (7 (v*7£  ) )  7? 


+  lower  order  terms 

0  T  0 

*  -(75  )  7y7£  +  transpose  +  lower  order  terms  • 

By  inserting  this  in  (5.18),  we  find  that  the  principal  part  of  the  constitutive  equation 
is  given  by 

(5.21)  [ft  +  U0 * V ] 2  "  g<x*t)7y(x,t)  , 

where 

t  35°  3cJ 

Mijmn  ’  _  _[  *l)kt 

Equations  like  (5.21)  can  also  be  derived  when  F  is  linearised  at  unsteady  flows  (§12). 


6.  SLIP  SURFACE  PROPAGATION  IN  PROBLEMS  PERTURBING  REST 

0  0 

In  problem*  which  perturb  rest  we  have  5  -  x,  x  »  0  and 

t 

(6.1)  T«  }  G(t  -  T)A1  [u(X,T)]dT  . 

A  Newtonian  part  of  the  stress  arises  when 

G( s )  •  yS{s)  +  g(s)  , 

where  6(s)  is  a  Dirac  measure  and  g(s)  is  a  smooth  kernel.  Thus 

t 

(6.2)  t  »  wA^ufx.t)]  +  j  g(t  -  T  )*1  [u(x,T))dx  . 


In  the  special  case  when  g(s)  is  in  exponential  form  (6.2)  reduces  to  the  model  (5.3)  of 
Jeffreys  with  u  -  hX^/X^.  He  introduce  u  here  and  elsewhere  only  to  notice  that  if  u 
is  small  the  underlying  dynamics  is  close  to  u  «  0  dynamics.  When  u  ”  0,  the  dynamics 
is  governed  by 

3t  t 

(6.3)  -  G(0)A1  [u(x,t)]  +  }  G’  (t  -  T)A1  [u(x,T)]dT  , 

and 

9u 

( 6. 3 >2  P  ^  +  Vp  -  div  T  -  0  , 

where 

(6.3) 3  div  u  -  0  . 

Equations  (6.3)  are  a  first  order  system,  linear  in  the  derivatives  of  p  and  the 
components  of  u  and  t.  If  G(t  -  T)  ■  const  exp{-(t  -  T)/X},  then  the  last  term  in 

(6.3)  ,  reduces  to  -t/X.  In  the  analysis  of  characteristics  using  the  Maxwell  model  the 
term  T/X  is  of  lower  order  and  it  does  not  enter  into  the  analysis  of  characteristics. 

In  general,  for  smooth  G' ,  in  (6.3)  the  integral  is  of  order  -1  in  t  and  +1  in 
x,  hence  of  order  0  as  an  operator  in  x  and  t,  and  is  thus  also  of  lower  order. 

This  shows  that  we  can  analyze  first  order  systems  for  general  kernels >  we  do  not  need 
special  models. 


R>>.  . . . 
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The  case  of  one-dimensional  shearing  motion  can  be  used  to  discuss  the  effects  of 
viscosity,  wave  speed,  wave  amplitude  and  classify  the  kernels.  In  this  case  (6.3)  reduces 
to 


3r_ 

at 


0(0)  +  /  G'(t  -  T) 


3x 


3u(  X,T  ) 
3x 


dx  , 


(6.4) 

We  may  write  (6.4)  as 

(6.5) 
where 


3u  3t 
0  3t  “  3x  • 


3g  3g 

3T+  §  37+  f-  -  °> 


1  -  “[ 


G(0) 

0 


1/pl 
0  - 


i  v  .* 

*v*v 

‘vV 

»,4  V  “ 


J  G'(t  -  T) 


3u(  X,T  ) 
3x 


dT 


The  eigenvalues 

(6.6) 


C  -  i  /G(0)/p 

of  B  are  the  wave  speeds  along  the  characteristics 

x  ±  ct  “  const. 

It  is  instructive  to  review  the  results  of  analysis  of  the  following  problem,  known  as 
Stokes'  or  Fayleighs'  problem,  for  a  viscoelastic  fluid.  A  fluid  occupying  a  half  space  is 
at  rest  for  t  <  0.  At  times  t  >  0,  the  boundary  of  the  fluid  at  x  -  0  is  made  to  move 
forward  with  a  constant  speed.  The  problem  may  be  described  by  the  following  equations 

i  +  *  '2 

3x2  —  3x‘ 

(6.7)  u(x,r)  -  0,  V  x  <  0  , 

u(x,t)  is  bounded  for  positive  values  of  x  and  t, 
r0,  t  <  0 


3u  32u  .  r 

p  JZ  *  11  —  +  '  g<t  ~  T) 


u(0,t)  -  H(t) 


1,  t  >  0  ’ 
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This  problem  has  been  studied  by  Narain  and  Joseph  [33]  when  u  *  0  and  u  ■  0  and 
by  Renardy  [37]  for  v  m  0  and  singular  kernels.  Other  authors.  Tanner  [46],  Straub  [45] 


This,  of  course,  corresponds  to  the  equation 

p _  pG1 (0) 

G(0>  Utt  '  G(0)2  Ufc  "  Uxx  ‘ 

Fran  this  one  expects  the  singularity  to  propagate  with  a  wave  speed  c  -  /g(0)/p  and  to 
decay  with  an  amplitude  factor  exp(tG' (0 )/2G(0 ) ) . 

Coleman  and  Gurtin  [6],  [7]  proved  that  if  the  acceleration  was  discontinuous,  then 
its  wave  speed  and  amplitude  must  be  given  by  these  expressions.  In  [33]  and  [37]  it  is 
shown  that  these  expressions  for  the  speed  and  amplitude  could  be  derived  for  the  problem 
(6.7)  without  assuming  discontinuities  of  acceleration.  The  demonstration  of  propagation 
in  [33]  applies  to  propagating  steps  in  displacement  (propagating  delta  functions)  as  well 
as  steps  in  velocity. 

If  a  Newtonian  term  is  included  in  the  constitutive  law,  G  contains  a  contribution 
|i£(0),  and  G(u)  at  m  behaves  like  0(1).  The  equation  becomes  parabolic,  and  smooth 
(analytic)  solutions  are  obtained.  As  u  +  0,  a  boundary  layer  forms  around  the  shock 
front.  This  was  shown  numerically  in  [46] .  An  analysis  of  this  boundary  layer  is  given  by 
Narain  and  Joseph  [33]. 

Renardy  [37]  studied  the  case  where  G  does  not  contain  a  S-function,  but  some 
weaker  singularity.  Specifically,  he  considers  the  kernels 


-GMt)  -  l 


-n*t  >  1 
,  «  >  J  ■ 


For  a  >  1,  G( 0 )  is  finite,  but  G'(0)  is  not.  The  asymptotic  behaviour  of  G(u>)  is 


given  by 


(6. 10) 


G(0)  ,  -2+1/a, 

B  +  0(»  >  • 


As  one  expects,  there  is  a  finite  wave  speed  c  -  /G(0)/p.  Solutions  are  zero  in  front  of 

00 

the  wave  and  nonzero  behind  it.  Across  the  wave,  however,  they  are  of  class  C  .  if 
a  <  1,  G(0)  is  infinite,  and  the  asymptotic  behaviour  of  G(w)  is  dominated  by  a  term 

1  /n-2 

0(u  '  )  or  Odn  itf/ai)  in  the  limiting  case  a  «  1.  The  wave  speed  is  infinite  and  the 

solution  becomes  analytic  everywhere  except  at  t  “  0. 

Similar  studies  can  be  done  for  small  perturbation  of  rigid  motions,  see  the  work  of 


Kazakia  [23]. 
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7.  SLIP  SURFACE  PROPAGATION  IN  NONLINEAR  SHEARING  PROBLEMS 


Coleman  and  Curtin  [6],  [7]  studied  simple  shear  flows  of  a  viscoelastic  liquid,  which 
involve  a  surface  across  which  the  acceleration  is  discontinuous*  They  showed  that  in  a 
material  with  a  smooth  memory  function  (no  Newtonian  viscosity),  these  discontinuities 
propagate  with  a  speed  c  _  /G(0)/p.  if  the  amplitude  is  smaller  than  a  critical 
amplitude,  it  decays  with  a  factor  exp( tG‘ (0 )/2G( 0) ) ,  if  the  amplitude  is  larger  than 
the  critical  amplitude,  it  will  reach  m  in  a  finite  time  ( "blow-up" ).  This  blow-up  of 
acceleration  waves  might  be  interpreted  as  development  of  a  slip  surface  for  the 
velocity.  No  such  result  has  bee  proven. 

Slemrod  [44],  Malek-Madani  and  Nohel  [28],  Gripenberg  [51]  and  Hattori  [16]  have 
studied  simple  nonlinear  models.  They  show  by  contradiction  using  Riemann  invariants  and 
the  method  of  characteristics  that,  for  suitably  chosen  initial  data,  a  global  C1 -solution 
for  the  equation  of  motion  cannot  exist.  This  may  mean  the  formation  of  a  slip  surface  for 
the  velocity.  This  is  suggested  in  particular  by  the  argument  of  Malek-Madani  and  Nohel 
[28],  but  no  full  proof  has  been  given.  Numerical  evidence  for  the  development  of  such 
discontinuities  was  found  by  Markowich  and  Renardy  [29]. 

Slemrod  [43],  following  a  suggestion  of  Coleman  and  Gurtin  [7],  applies  his  result  to 
give  an  explanation  of  some  of  the  various  instabilities  of  shear  flows  collectively  called 
melt  fracture.  It  is  necessary  that  we  interpret  the  proved  loss  of  regularity  as  the 
formation  of  a  slip  surface  for  the  velocity.  Melt  fracture  is  an  instability  of  flow  of 
molten  polymers  or  polymeric  solutions  down  capillaries.  In  the  experiments  (Tordella 
[47]),  the  polymer  is  forced  down  the  capillary  by  high  pressure.  Extrudates  leaving  the 
capillary  which  at  lower  shear  rates  are  smooth  and  continuous,  become  rough  (shark  skin 
effect),  irregular  and  ultimately  disintegrate.  There  are  different  explanations  of  the 
different  types  of  instability  which  can  occur.  These  are  reviewed  in  the  paper  of  Petrie 
and  Denn  [35].  None  of  the  explanations  can  be  regarded  as  established.  The  mechanism 
proposed  by  Slemrod  has  some  possibilities  in  the  case  where  fracture  is  associated  with  a 
stick-slip  phenomenon.  There  is  some  controversy  about  the  presence  or  absence  of  slip  in 
experiments.  If  it  does  occur  periodically,  as  it  might  in  Slemrod' s  theory,  it  would  be  a 
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candidate  for  the  explanation  of  the  wavy  surfaces  shown  in  the  pictures  of  Tordella 
[47].  It  should  be  noted  that  the  theory  of  Hunter  and  Slemrod  [17]  requires  an  entirely 
different  type  of  shear  stress  to  explain  the  type  of  hysteretic  melt  fracture  which 
Tordella  calls  ripple. 

It  would  be  interesting  to  have  the  conditions  under  which  slip  surfaces  for  the 
velocity  might  develop  from  weaker  slip  surfaces  or  from  smooth  data. 


•»  V 
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8.  CLASSIFICATION  OF  EQUATIONS  FOR  FLOWS  PERTURBING  UNIFORM  MOTION 

Ultman  and  Dann  [49]  consider  the  equations  for  two-dimensional  steady  flow  of  an 
upper  convected  Maxwell  fluid.  They  linearize  at  a  motion  with  uniform  velocity  and  zero 
stress,  and  they  show  that  these  linearized  equations  change  type  when  a  viscoelastic 
"Mach"  number 


U 


*> 

i! 

V 


(8.1)  M  -  - 

c 

exceeds  one.  Here  U  is  the  velocity  of  the  unperturbed  uniform  flow,  and  c  is  the  wave 
speed  for  propagation  of  singularities  as  considered  in  $6:  c  “  /G(0)/p  »  /n/(pX). 

Ultman  and  Denn  attempted  to  correlate  some  experimental  observations  of  D.  F. James 
[18]  with  this  change  of  type.  James  observes  a  sudden  change  in  the  slope  of  the  heat 
transfer  curve  as  a  function  of  velocity.  This  happens  at  a  critical  velocity,  which,  for 
the  polyox  solution  used  by  James,  was  about  1  cm/sec.  It  is  not  clear  from  the  graphs  how 
abrupt  this  change  in  slope  is,  but  there  is  a  change  of  slope.  Ultman  and  Denn  also 
suggest  that  the  transition  from  subcritical  to  supercritical  flow  might  explain  abrupt 
changes  in  the  drag  coefficient  observed  by  A.  Fabula  [13].  Again,  the  idea  is  that  the 
critical  velocity  at  transition  is  the  wave  speed  c.  They  make  an  estimate  of  c  from  a 
molecular  theory  and  correlate  this  prediction  with  the  data  of  James.  Of  course  any  such 
estimate  can  at  best  be  expected  to  give  an  order  of  magnitude,  since  the  fluids  used  in 
experiments  are  not  really  Maxwell  fluids. 

M.  Luskin  [27]  studies  the  equations  of  Ultman  and  Denn  as  a  first  order  system.  He 
reduces  them  to  canonical  form  and  investigates  characteristics.  The  streamline  is  a 
double  characteristic,  two  characteristics  are  always  complex,  and  the  remaining  two  are 
complex  for  M  <  1  and  real  for  M  >  1.  Each  characteristic  has  associated  with  it  a 
canonical  variable,  which  is  a  linear  combination  of  the  two  velocity  components,  the  three 
components  of  the  extra  stress  and  the  pressure.  This  means  that  two  of  these  variables 
can  be  discontinuous  across  streamlines,  and  two  others  can  be  discontinuous  only  for 
M  >  1. 

We  shall  now  show  how  results  similar  to  those  of  Ultman  and  Denn  apply  to  general 
constitutive  models  without  "viscosity".  In  linearization  at  uniform  motion,  the  extra 
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stress  is  given  by 


T  -  /  G(t  -  T  )J[u(£,T)]dT  , 


X,  -  U(t  -  T) 

5  -  L1  ) 


He  assume  G  is  smooth,  positive  and  monotone  decreasing.  By  differentiating  (8.2)  with 
respect  to  t,  holding  x  fixed,  we  find 


3t  3t  t 

+  U  3^- -  G(0)A[u(x,t)]  +  /  G*  (t  -  T)A[u(£,T)]dT, 


where  the  last  term  is  of  lower  differential  order.  The  leading  term  is  the  same  as  for 
the  Maxwell  model,  if  G(0)  is  replaced  by  n/A. 

It  can  be  shown  that  the  change  of  type  is  primarily  associated  with  the  behaviour  of 
the  vorticity.  Since  we  study  motions  in  the  plane,  the  vorticity  curl  u  has  only  one 
component,  which  we  denote  by  a.  He  take  the  curl  of  the  equation  of  motion,  apply  the 
operation  curl  div  to  (8.3)  and  combine  the  two.  In  steady  flow  this  leads  to 

.2  .2  t 

(8.4)  (M2  -  1)  - /  G’(t  -  T)Aa(£)dT  . 

3x^  3Xj  — 

The  right  side  is  again  of  lower  order.  The  left  hand  side  is  elliptic  when  M  <  1  and 
hyperbolic  when  M  >  1.  The  elliptic  roots  found  by  Luskin  [27]  correspond  to  the  elliptic 
equations  a  «  -&i>  expressing  the  vorticity  in  terms  of  the  stream  function. 
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The  flow  of  an  upper  convected  Maxwell  fluid  Is  governed  by  the  following  system  of 
equations 

du 

p  —  +  7p  -  dtv  T  ■  0,  div  u  «  0  , 

<9‘1)  at  _ 

T  +  x[“  -  Vu  T  -  T  VuTj  ■  nA  “  pXA  . 

We  consider  two-dimensional  flows.  In  view  of  the  applications  in  Chapter  10,  we  want  to 
write  the  equations  in  both  Cartesian  and  polar  coordinates.  In  cartesian  coordinates,  we 


r  0 

U  -  (u,V),  T  «  lT 


We  then  obtain  the  equations 


a  +  up  +  vo  -  2ou  -  2tu  -  2pu  -  -o/X  , 

t  x  y  x  y  x 

t  +UT  +  VT  -YU  -0V  -  UtU  +  V  )  •»  -T  A  t 

t  x  y  y  x  y  x 

y.  +  uY  +  VY  -  2tv  -  2yv  -  2uv  ”  -Y/X  , 

't  x  y  x  y  y 


p(u  +uu  +  vu)+p  -  o  -  T  “0 
t  x  y  x  x  y 

p(vt  +  uvx  +  wy)  +  py  -  Tx  -  Yy  -  o 
ux  +  vy  -  0  . 


The  subscripts  denote  differentiation.  In  polar  coordinates,  we  denote  by  u  and  v  the 
radial  and  azimuthal  velocities  and  by  o,t,y  the  components  of  the  stress  in  polar 
coordinates.  We  then  obtain  the  following  system 


vo0  u0  9 

a  +  uo  - - 20 u  -  2t - 2yu  -  -  t-  , 

t  rr  r  r  r  X 


vTe  Yue 

T  +  UT  +  -  -  -  -  <3V 

t  r  r  r  r 


-  „(v  +-A).-i.ax.iu, 

v  r  r  '  X  r  r 


\  +  uYr  +  ~T  -  2TVr  '  ^  IT  ~  2V  7~  “  “  X  "  T1  +  ?  (Y  +  U  »  ' 
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plu.  +  uu  +  v  — I  +  D  -  a  -  —  -  p  —  +  - 

v  t  r  r  y  *r  r  r  r  r 


vv  p8 


/  _ 0  0 

P  [  V.  +  UV  +  - 1  +  -  -  T  -  -  ! 

vt  r  r  '  r  rr 


puv  2JL 

r  r 


8  u 

u  +  —  * - 

r  r  r 


The  terms  on  the  right  of  (9.3)  and  (9.4)  are  of  lower  order  and  do  not  enter  into  the 

analysis  of  characteristics*  The  systems  on  the  left  are  identical  if  we  make  the 

3  3  3  1  3 

identifications  t—  -  t-~  and  v—  ~  —  -rr. 

3x  3r  3y  r  30 

He  consider  steady  solutions  of  (9.3)  and  (9.4).  Thus  we  put  the  time  derivatives 
equal  to  zero  and  introduce  3  -  (u.v.ct ,t ,y ,p) .  Equations  (9.3)  and  (9.4)  are  of  the 


respective  forms 


SSx  +  S3y 


A  A  2ft  A 

(9.6)  Ag^  +  B  —  »  f  . 

Characteristics  of  (9.5)  are  determined  by  the  equations 

(9.7)  det(~£  A(g)  -  B(g) )  =  0  . 

With  a  »  -  — this  leads  to 
dx 

2  2  2  2  2 

(9.8)  (au  +  v)  (1  +  a  )(2ax  -  p(au  +  v)  +ao+Y+(1+a)p]«0, 

Hence  the  streamlines 

(9.9)  ^  u  -  v  -  0 

dx 

are  double  characteristics.  There  are  two  imaginary  roots  a  =  +  i.  Finally,  the  last 
bracket  yields 


u  +  a  -  pu 


(U  +  Y)(u  +  o)  +  pv2(p  +  0)  +  pu2(p  +  y  )  1  V-> 


(u  +  0  -  pu  ) 


These  characteristics  are  real  or  complex  depending  on  whether  the  argument  of  the  square 
root  is  positive  or  negative.  ThiB  argument  of  course  depends  on  the  solution.  It  should 
be  kept  in  mind  here  that  the  integral  form  of  the  Maxwell  model  imposes  certain 


.v.’.  ,%  ,>  .V  .%  „\  .1  s 


restrictions  on  ths  stresses  (see  Chapter  5).  These  roots  are  the  ones  which  exhibit  a 
change  of  type. 


The  imaginary  roots  ^  ■  +  i  arise  from  the  identity 


(9.11) 


?  -  -V**  , 


where  C  is  the  vorticity  and  *  is  the  steam  function.  It  is  easy  to  show  that 

2  S2C  2  32r  j2r 

(9.12)  (U  +  o  -  DM  )  — |  +  (U  +  Y  -  pv  )  — |  +  2(t  -  puv)  3x^y  -  lower  order. 


This  shows  that  the  interesting  characteristic  roots  (9.10)  which  can  change  type  are 
associated  with  the  vorticity.  Using  (9.11),  we  find  that  the  quantity  on  left  of  (9.12) 
with  ♦  replacing  C,  is  harmonic  to  within  terms  of  lower  order.  In  polar  coordinates 
the  equation  for  chgaracteristics  becomes 


(9.13) 


r  ■  r.h.s.  of  (9.10)  , 


with  u,v,T,Y,a  Interpreted  as  the  physical  components  of  velocity  and  stress  in  polar 
coordinates . 


■'.N'VWV 
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10.  CHANGE  OF  TYPE  IN  SHEAR  FLOW,  EXTENSIONS  FLOW,  SINK  FLOW  AND 
CIRCULAR  COUETTE  FLOW  OF  AN  UPPER  CONVECTED  MAXWELL  FLUID 

In  this  section  we  consider  the  problem  of  hyperbolicity  of  the  linear  equations 
perturbing  some  special  solutions  of  (9.1).  The  characteristic  equations  (9.10)  and  (9.13) 
are  useful  for  this.  To  compute  the  characteristics  of  the  perturbation  we  need  only  to 
replace  the  u,v,t,<j,y  in  these  formulas  with  the  special  values  that  (9.1)  requires  for 
special  flows. 

(i)  First,  we  again  consider  the  uniform  flows  discussed  in  <8;  u  *  U, 
v»t»Y«o«0.  For  these  we  find  from  (9.10)  that 


where  »  U^/C^,  C^  ”  n/Xp(«  ' )  is  the  wave  spread.  The  characteristics  are 

straight  lines  which  start  as  lines  perpendicular  to  the  flow  at  M  »  1  and  tilt  more 

toward  the  free  stream  as  M  >  1  is  increased. 

(ii)  Our  second  application  is  to  shear  flow.  We  find,  using  (9.3),  that 

2 

(10.2)  u  “  <y,  t  **  n*>  a  »  2nXic  ,  y  =  v  »  o  . 

where  the  shear  rate  <  is  a  constant.  Inserting  the  fields  (10.2)  into  (9.10),  we  find 
that 


(10.3) 


dx 


The  vorticity  is  hyperbolic  outside  a  strip  defined  by 


\  f 


(iii)  In  extensional  flow,  we  find,  using  <9. 3),  that 


(10.5) 


u  “  ax,  v  =*  -sy,  x  -  0 ,  a  *  2ns/(1  -  2Xs),  Y  m  -2ns/(1  ♦  2Xs)  , 


where  0  <  s  <  1/2X.  The  stretch  rate  a  is  a  positive  constant,  small  enough  to  keep  o 
bounded  and  positive.  The  unbounded  0  at  s  »  1/2X  is  one  of  the  many  undesirable 
properties  the  upper~convected  Maxwell  model.  Inserting  the  field  (10.5)  into  (9.10),  we 
find  that 


(10.6) 


dx  U  2  2 

7^x7"  ps  x 


- -  1  2  2 

.  .,2  2  a  b 

1  -  4X  s 


2  2 

{*_  +  iL.  .  1}1-/2  , 


ps  (1  -  2Xs) 


-,  b2  - 


ps  (1  +  2Xs) 


We  get  real  characteristics,  obeying  (10.6),  outside  the  ellipse 

2  2 

VV1- 

a  b 

(iv)  Our  fourth  application  is  to  sink  flow.  Here  we  must  work  with  the  polar 
equations  (9.4)  and  (9.13),  For  sink  flow 


(10.7) 


u  -  -Q/r,  Q  >  0  is  the  sink  strength 


v  -  x  -  0  . 


From  (9.4),  we  find  that 


We  find  that 


(  10.8) 


From  (9.4)3,  we  f*n<J  that 


5®  +  2£_£5+2E,0. 

dr  r  XQ  r 


o  »  2n  Q/r 


51  .  21  .  SL  -  SH.  -  0 

dr  r  XQ  r 


We  next  Introduce  ♦  through  the  change  of  variables  Y  “  ~U  *  V$AQ,  where 


SSI 


5i  _  *[2  +  r  1  +  r 

dr  ,Lr  XQJ 


We  require  that  $  be  bounded  as  r  ♦  0  and  find 


(10.9) 


dr  , 


♦  -  -rV  /2X°  /  — — 

r 

where  the  last  factor  is  an  indefinite  integral.  For  small  r,  we  find  that 

2 

♦  '■  -r  log  r  >  0  . 

Obviously  4  is  poaitive  at  every  stationary  point.  The  only  stationary  point  is  at 
r  »  •  and  4  increases  monotonies lly  from  zero  to  XQ. 

Turning  next  to  (9.13),  using  (10.7-9),  we  find  that  real  characteristics  exist  and 
are  given  by  (9.13)  when 

no.io)  (v  +  y)[£S —  in  +  «)]■§  [— — ~  2n^  -  u]  >  o  . 

r  “  r 


The  condition  for  hyperbolicity  is  satisfied  if 

PQ  >  2n  » 

that  is,  when  the  flow  rate  is  large.  In  this  case  the  flow  is  hyperbolic  in  the  circle 

(10.11)  r  <  J  e£-  ~  . 

The  differential  equations  for  the  net  of  characteristics  covering  the  hyperbolic 
circle  at  the  origin  la 


r  4®  _  +  !  H  +  *  -  +  /  ..  -U»/x2  -  . 

^  /  Pu2-(M0|  /  p,g2_r  *12  _  „ 

2  w 
r 

(v)  Our  fifth  application  is  to  Couette  flow  outside  a  rotating  circular  cylinder  of 
radius  a.  We  suppose  that  the  fluid  sticks  to  the  rod,  and  the  rod  rotates  with  an 
angular  frequency  SI.  He  assume  that  v  “  v(r),  u  ■  a  ■  0.  Then  using  (9.4)5  and  (9.4)2 
we  find  that 

2 

v  »  fla  /r  , 

T  -  -2a2Wr2  , 


and  from  (9.4 )3 


r  -  2t2/u  . 
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Fran  (9.13),  we  find  that 


To  have  real  characteristics  it  is  necessary  and  sufficient 


2,  def  a4!)2 

E(r  ’  m  -22 
c  r 


a"fl~  4X2c2-i  ,  v  . 

-2-2(1-  —r~)  -  1  >  0 


Moreover,  Z(m) 


dE _  a4!)2  r8X2c2  .1 

2  2  4  l  2  “  ’J 

dr  c  r  r 

2  2  2 

first  increases,  then  decreases,  with  a  single  maximum  at  r  ■  8X  c  given  by 

1  , 


T.  -  E(8X2c2) 


4n2 
a  ii 

.,,2  4 

16X  c 


which  is  positive  when 

(10.12) 

If  (10.12)  holds  then  there  are  two  values 


a4 


8X2c2 


.,/t: 


16X  c 

o2  4 
fl  a 


at  which  S(r^)  -  0.  if 
c 


a2  >  4X2c2  +  %  , 


then  E(a)  >  0  and  the  vorticity  is  hyperbolic  with  real  characteristics  given  by  (9.13) 
in  the  annulus 


2 

a 


< 


8X2c2 


1  -  J  1  - 


16X2c4 

n2a4 


2 

a 


<  4X2c2  +  ^  , 

n 


If  (10.12)  holds  and 


then  £( a)  <  0  and  the  vorticity  is  hyperbolic  with  real  characteristics  given  by  (9.13) 
in  the  annulus 


(vi)  M.  w.  Johnson  (private  communication)  has  studied  the  flow  of  an  upper  convected 
Maxwell  fluid  between  eccentric  rotating  cylinders.  He  assumes  a  small  gap  and  uses  the 
lubrication  approximation.  His  analysis  shows  that  change  of  type  occurs  at  high 
Weissenberg  numbers. 


11.  CLASSIFICATION  OF  THE  QUASILINEAR  SYSTEM  OF  EQUATIONS 
GOVERNING  STEADY  FLOWS  OF  A  CLASS  OF  OLDROYD  MODELS 

The  Oldroyd  models  under  discussion  are  given  by  (5.4)  and  (5.5).  When  a  «  1, 
model  reduces  to  the  upper  convected  Maxwell  model  which  we  studied  in  $9  and  10. 

We  have 


This  gives  rise  to  the  following  quasilinear  system 

ucx  +  vay  +  t(vx  -  uy)  -  »t2oux  +  T(Uy  +  vx)J  -  2uux  -  -a/\  , 

urx  +  vry  +  1  (o  -  Y)(uy  -  vx)  -  |  (a  +  Y)(uy  +  vx)  -  u(uy  +  vx)  -  -t A  , 

(11.1)  uY  +  VY  +  T(U  -  v  )  -  *(2yv  +  T(u  +  v  ))  -  2yv  -  -  *  , 

*  y  y  x  y  y  *  yA 

p(uux  +  vuy)  ♦  px  -  Ox  -  Ty  -  0  . 

p(uvx  +  wy)  +  Py  -  Tx  -  Yy  -  o  , 

u  +  v  »  0  . 

X  y 

The  analysis  of  characteristics  follows  the  one  used  in  (9  exactly.  We  find 
characteristic  directions,  a  -  -dy/dx  from  the  equation 

(11.2)  (1  +  <*2)(<*u  +  v)2{p(au  +  v)2  +  (^2-)(a2  -  1)  -  2to  -  (a2  +  1)(p  +  a(^p))  -  0 
corresponding  to  (9.8).  We  can  prove  that  the  vorticity  changes  type  from  elliptic  to 
hyperbolic  when  the  sign  of 

(11.3)  (Pu2  +  ^  (1  -  a)  -  ~  (1  +  a)  -  Ml  1(1  +  a)  ^  +  (a  -  1)  |  +  u  “  Pv2]  +  (puv  -  t) 


changes  from  negative  to  positive 


( 11.6)  -pu  +  u  +  “  (  i  +  a)  +  T  (a  -  DC  +  2(-puv  +  T  )c 

'  2  2  xx  xy 


f-pv2  +  u  +  7  (a  -  1)  +  £  (a  +  1)]c  +  Jt.o.  -  AC  +  2BC  +  CC 

1  2  2  J  w  def  xx  xy  yy 


This  a  mixed  elliptic-hyperbolic  second  order  equation.  It  is  hyperbolic  for 


A  «  B  -AC  >  0,  l.e.  exactly  when  (11.3)  holds  (see  the  previous  discussion). 


nation  for  vortlclty,  unsteady  case. 


In  unsteady  planar  flow  the  equation  for  vorticity  is 


2 

(11.7)  P(f£  ♦  <u-V)C)  +  X[p  +  2p(u»V)  || 


+  (pu2  -  u  -  f-  (1  +  a)  +  ^  (1  -  a))  ^  + 


2  +  2<PUV  “  T) 


(pv2  -  m  +  f(1  -  a)  -  ?•  (a  +  1))  -  X  [i.ol  . 


The  term  in  brackets  on  the  right  hand  side  is  of  lower  order  and  does  not  involve  any  time 


derivative.  The  first  term  on  the  left-hand  side  is  irrelevant  for  the  analysis  of  the 


type  of  equation  (11.7),  but  it  shows  how  (11.7)  reduces  to  the  unsteady  Newtonian  case 


whan  X  «  0  (recall  that  u  »  •j-). 


Equation  (11.7)  is  of  evolution  type  (hyperbolic)  when 


(11.8)  X2T2  -  [n  -  x(-j  (1  -  a)  -  ~  (1  +  a) ) ]•  [n  -  X  (-|  (1  -  a)  -  ^  (1  +  a))]  <  0  , 


(11.9) 


X[j  (1  -  a)  -  ~  (1  +  a)]  -  n  <  0 


i  s 


12.  QUASILINEAR  SYSTEMS  FOR  SIMPLE  FLUIDS  WITH  FADING  MEMORY  OF  THE  COLEMAN-NOLL  TYPE 

The  determinate  stress  T  in  a  simple  fluid  is  given  by  an  isotropic  functional  of 

T 

the  history  of  the  relative  Cauchy  strain  G(s)  -  F  (t  -  s)F  (t  -  s)  -  1,  i.ee 

*  *t  *t  ■ 


M2.1)  qtqt  -  p[  g  G(s)gTj 

s»0~ 

for  all  constant  orthogonal  tensors  g.  By  taking  the  material  derivative  of 
obtain 


(12.2) 


(12.1),  we 


where  F^  denotes  the  first  functional  derivative  of  P.  We  have  alreadynoted  in  §5  that 
different  choices  for  the  domain  of  the  linear  functional  P^ [QGQ  | • ]  lead  to  different 

representations  of  P^ .  If,  following  Coleman  and  Noll,  we  assume  that  the  functional  is 
2 

defined  on  a  weighted  L  -space,  we  obtain  an  integral 

dG  «  dG(s) 

O2.3)  pjG  I  ^  ]  -  J  K(s,G)  ds  . 

0  « 


Here  K(a,G)  is  a  fourth  order  tensor  depending  on  s  and  on  the  value 
{G(o),  0  <  a  <  «•}.  For  the  following,  we  assume  that  X  and  its  firBt  derivative  (as 
functions  of  s)  are  integrable.  In  particular,  it  follows  that  K  is  uniformly  bounded 
in  a. 

The  isotropy  condition  can  be  written  as 

(12*4)  >W'S>  -  CaiCbjPck°dlKabcd(8'22ST)  • 

This  consequence  of  isotropy  does  not  appear  to  lead  easily  to  very  explicit 
representations  for  K  in  situations  where  G  is  not  confined  to  special  motions.  Of 

m 

course  X  is  syometric  in  the  first  two  indices,  and  only  the  symmetric  part  in  the  last 
two  indices  enters  into  (12.3). 

We  next  note  that  the  material  derivative  of  G  is  given  by 


II 


\ 

1 

1 

h 
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(12.5) 


■=-  G(S)  -  —  (F  ( t)F  (t  -  s)F(t  -  8  )  F  ( t  >  > 

dt  “  dt  «  -  *  » 


-  -F'T(t)FT(t)G(8)  “  G(s)F(t)F_1 (t)  -  G(s) 

-  “  =  “  ”  ds  “ 


T  dS<®> 

-L  G(s)  -  G( a )L  -  — —  , 


where  J(x,t)  »  Vg(x,t)  is  the  present  value  of  the  velocity  gradient.  Hence,  we  find 

(12.6)  /  Klj)tl<.,G)  d 8  -  -  J  (Kijkl  ♦  Kijlk)Gpl(B)ds  •  Lpk(t> 

-  dg(s) 

"  I  W'S'  —  d8  • 


It  follows  that 


(12.7) 


"dt  “  “ijkp'S'^pk1*'0  +Nij(“'t)  ' 


(12.8) 


•W*0  -  •{  <Kijki +  Kijik,Gpi(e)d8 ' 


“  *5  <s)  "  a 

(12.9)  «lj(*.t)  -  /  Kljkl("'S(S,)  Si -  d8  -  "J  Gkl(8)  d?  Kijkl(B.G<B)’dB  * 

■Hie  coefficients  Mijkp  ot  Lpk  and  the  terms  (  jj)  are  of  lower  order  and  j <  S't ) 
does  not  enter  into  analysis  of  problems  of  change  of  type. 

He  may  write  (12.9)  as 


(12.10) 


“  8m  D,  +  A,....  ft.  +  N  , 
dt  ijkp  kp  ijkp  kp  ij 


Mijkp  “  si jkp  +  Aijkp  ' 

and  8  is  symmetric  and  A  skew  symmetric  in  kp.  D  and  Q  denote  the  symmetric  and 
skew  symmetric  part  of 


In  its  general  form,  (12.10)  and  the  equations  of  motion  are  a  auacilinear  system  of 


first  order  in  the  derivatives  of  p,  y,  £.  We  could  write  out  conditions  for 
evolutionarity  or  for  change  of  type  in  steady  flow  as  before.  However,  we  can  not  in 
general  isolate  an  equation  for  the  vorticity,  as  we  did  for  the  Oldroyd  models. 

He  can  identify  a  class  of  models  which  is  more  special  than  the  completely  general 
equation  (12.10),  but  much  more  general  than  the  three-constant  Oldroyd  model.  Let  us 
assume  that 


(12.11) 


;  «  1  ({  p  +  6  P  +  6  p  +  {  p  ) 

ijkp  2  ik  jp  jk  ip  ip  jk  jp  ik* 


where  £  is  any  second  order  tensor,  expressible  by  integrals  of  the  type  (12.8).  He  need 
no  assumptions  at  all  on  the  anti-symmetric  part  ^i j^p*  Using  (12.11)  we  may  reduce 
(12.10)  to 


PikDkj  +  PikDki  +  Aijkpnkp  +  Nij 


(12.12) 


a“  T  1  T 

—  -  PD  +  DP  +  -  (A*G  +  (AG)  )  +  N 
dt  "  2  -  *  “  - 

where  jkp  i8  symmetric  in  i j  and  skew  symmetric  in  kp . 

The  three  constant  Oldroyd  model  arises  from  (12.12)  with  special  choices  for  the 
tensors  $,£  and  Hi 


(12.13) 


-2Aijkl  "  aikTJl  +  SjkTil  ”  SilTkj  "  SjlTki  ' 


ik 


•Tik  +  X  sik  ' 


ij 


7  Tij  • 


He  now  demonstrate  that  the  quasilinear  system  associated  is  (12.12)  is  expressible  in 


terms  of  the  vorticity  5  “  curl  u.  The  equations  of  motion  are 
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(12.14) 


du 

p  ~  "  -Vp  +  div  T ,  div  u  »  0 
at  * 


He  apply  the  operations  curl  and  —  to  (12.14),  and  find  that  to  leading  order 

dt 


(12.15) 


d2  d- 

p  ' — —  C  ■  curl  div  —  +  . . .  . 

dt2  “  dt 


I 

ti 


I 


The  question  is  now  whether  the  right  side  of  (12.15)  is  expressible  in  terms  of  t ,  to 
leading  order.  Clearly  fl  is  expressible  in  terms  of  ;<  so  we  need  no  restrictions  on 
A.  Then,  working  the  part  with  tensor  P,  we  calculate 

'abiVj^kj  +  P3k°ki)  “  W'ikVAj  +  Vb3jDki]  +  ••• 


33u, 


1  e  p  v2_\+1e  - * _ 

2  abi  ik  2  abi  jk  3x^3x^3x^ 


+  . . 


a 


Noting  now  that 


V  u  “  -curl  C  , 


V*( 


we  get 


~  2  EabiPik3b(CUrl  ^}k  +  2  Pjk3j3kCa  +  **• 


It  follows  now  that  to  leading  order  (12.15)  is  a  second  order  quasilinear  system  of 
equations  for  the  components  of  the  vortlcity. 
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